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In this short letter, we report on lepton flavor violating Higgs decay into [if + r/2 
under the scheme of a generic supersymmetric standard model without R-parity and 
list important combinations of R-parity violating parameters which can give signifi- 
cant branching ratios in experimental accessible region. We impose other known ex- 
perimental constraints on the model parameters, considering both with and without 
a sub-eV neutrino masses bound on mass contributions from the R-parity violating 
parameters. In our analysis, the branching ratio of h° — > \xf + Tp can exceed 1CP 5 
within admissible parameter space, while 10 -2 is possible when the sub-eV neutrino 
masses bound is lifted. 
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I. INTRODUCTION 

In standard model (SM), lepton number of each flavor is separately conserved. However, 
it is well known now from the observation of neutrino oscillation experiments that lepton 
flavor conservation should be violated which also gives hints of new physics beyond 
the standard model. Therefore, various lepton flavor violating (LFV) processes have gotton 
much attention and been analyzed. For example, in the framework of supersymmetry, which 
is no doubt the most popular candidate theory for physics beyond SM, flavor violating 
processes such as r — > /ry, r — > fiX, r — )■ /it/ and r — > fififi etc. are discussed in minimal 
supersymmetric standard model (MSSM) [2] and other supersymmetric extensions jsl, Q]. 
Specifically, for LFV Higgs to /if + r/2 decay which we put our focus on in this letter, has 
also been studied generally in the general two Higgs doublet model jlj], which allow lepton 
flavor violation, and in the MSSM with right-handed neutrinos fjj]. Though supersymmetry 
is so far in 



pointed out 



ack of experimental evidence and hence faces stringent challenges, it is 
8| that there are still rooms for (minimal) supersymmetric standard model to 
accommodate existing experimental results. For example, large mass spectrum around or 
beyond ITeV for majority of supersymemtric particles is possible; by allowing non-universal 
soft breaking masses or more free parameters can also relax the constraint from experiments. 

In many studies of the supersymmetric standard model, a discrete symmetry called "Im- 
parity", which keeps baryon and lepton number conservation, is often imposed by hand to 
prevent proton decay and make the lightest supersymmetric particle (LSP) a dark matter 
candidate. From the theoretical point of view, R-parity is ad hoc, not motivated so long 



Oth- 



as the phenomenological (minimal) supersymmetric standard model is concerned 
erwise, a generic version of (minimal) supersymmetric standard model without R-parity, 
which is our focus here, has the advantage of a richer phenomenology and giving neutrino 
masses and mixings without the need to introduce any extra superfield. Beyond neutrino 
masses, the key part of R-parity violation of interest is lepton flavor violation, because it 
can be observed in processes involving only SM particles while SM itself does not give rise 
to this kind of processes. Higgs to /if + r/2 decay is our example at hand. In the framework 
of R-parity violating (RPV) supersymmetry, there have been some studies 1CH12| on the 
issue of lepton flavor violation. However, such studies either were limited to particular types 
of R-parity violation or did not study Higgs to /if + r/2. 
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Recently both ATLAS and CMS reported discovery of a boson state with mass near 125- 
126 GeV jisifl^l while the Tevatron CDF and DO collaborations are observing a statistically 
significant excess in their bb + W or bb + Z channels [jjj]. The boson is essentially compatible 
with a SM-like Higgs, while more data is needed to fully pint down its nature. The flavor 
violating decay h° — > \ir + r/2 is of particular interest at this moment. 

In this short letter, we present the first results of a comprehensive study on the generic su- 
persymmetric standard model, without R-parity, highlighting; cases of most interest. Details 
of our studies will be reported in another publication [16 1. 

In following section, we summarize the basic formulation and parameterization of the 
generic supersymmetric standard model. The neutral scalar mass matrix is particularly 
discussed, including the key quantum corrections up to two-loop order. Then we sketch our 
calculations and present our numerical results in section 3. Note that during our analysis we 
made no assumptions on RPV parameters hence not losing general validity while keeping 
Higgs and the mass spectrum of the other particles within experimental constraints. Finally, 
we conclude this letter and give an account of upcoming works in the last section. 



II. GENERIC SUPERSYMMETRIC STANDARD MODEL (WITHOUT 
R-PARITY) AND HIGGS MASS MATRIX 



With the superfield content of the MSSM, the most general renormalizable superpotential 
can be written as 



W = e ab 



where (a, b) are SU(2) indices with e i2 = — e 2 i = 1, (i,j,k) are the usual family (flavor) 
indices, and (a, (3) are extended flavor indices going from to 3. Note that A is antisymmetric 
in the first two indices as required by SU(2) product rules while A" is antisymmetric in the 
last two indices by SU(3) C . The soft SUSY breaking terms, which can be written as it 
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follows: 



soft 



: £abB a H°;L h n + e 



~ab D OL — u — 
+ Cab 



aft 



AmH b u Uf + A?H%Q<>Df + 



+ h.c. 



M] - - M 2 - - M 3 

+ + -fWW + + h.c. , 

where L*m~L is given by a 4 x 4 matrix with zeroth components, m 2 



(2) 



Lqq corresponds to fh Hd 
in MSSM while m| ofc 's give new mass mixings. 

The above, together with the standard (gauged) kinetic terms, describe the full La- 
grangian of the model. We have four leptonic superfields L, which contain the components 
of fermion doublet as 1° and l~, while their scalar partners as 1° and /~. For convenience, we 
choose a flavor basis such that only Lq bears a nonzero vacuum expectation value (VEV) 
and thus can be identified as the usual Hd in the MSSM. However, one should keep in mind 
that Hd now may contain partly the charged lepton states. The parameterization has the 
advantage that tree level RPV contributions to fermion mass matrices are described com- 



pletely by the /ij, Bi, and rn^- parameters, which are wef 



just very conservative neutrino mass bounds imposed 17 



constrained to be small even with 



181 ] . The tree level mass matrices 



for the scalar fields also have minimal extra terms from R-parity violation [17] . 

For the mass matrices for the (colorless) scalars, we have five neutral complex fields, 
from H u and four L a superfields, and eight charged fields. Explicitly, for the neutral part 
we write the (1 + 4) complex fields in terms of their scalar and pseudoscalar parts, in the 
order , Zq, ^Dj t° form a full 10 x 10 (real and symmetric) mass-squared matrix. 

As for charged scalars, the basis {h^ Jo Jt Jt } * s use d to write the 8x8 

mass-squared matrix. The exact form of these tree level matrix elements can be found in 



17|. 



With all the RPV terms, the physical scalar states are now mixture of Higgses and 
sleptons. The RPV terms provide new contribution to scalar mass matrices and hence 
Higgs mass. In addition, third generation quarks and squarks could play an important role 
in radiative corrections to the Higgs sector and hence should be included. Accordingly, we 
implement complete one-loop corrections, from Ref . 19) , to matrix elements directly relating 
to Higgses (CP-even, CP-odd and charged ones as well) during our computation. Moreover, 
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light Higgs mass issue is quite sensitive and need to be delicately treated, therefore we 



include an estimation 20] of the key two-loop corrections in light Higgs related elements 
also 1 . Note that radiative RPV corrections are typically too small to be taken into account, 
thus we study tree level RPV effects only. 

III. THE CALCULATIONS AND NUMERICAL RESULTS 

In this section, we describe most important RPV parameter combinations which are 
expected to have largest branching ratios without violating existing experimental results. 
There are several different sources which can give constraints on our RPV parameter setting. 



Among which, the one from indirect evidence of neutrino mass, i.e. £V m Ui < leV [21] is 
quite crucial. However, since it has not been completely ruled out for neutrinos having mass 
larger than leV, we give some comments on branching ratios as reference under the condition 
that neutrino mass constrained only by the solid bounds, i.e. m Ue < 3eV, m Ufl <190keV and 
m Ur < 18.2MeV as well. 

In our numerical computation, we deal directly with mass eigenstates. We put all the tree 
level mass matrices into the program. In the case of Higgs masses, we include as well the 
one-loop and two-loop corrections to relevant mass matrix elements as mentioned in section 
2. The mass of the Higgses (and other sparticles) needed in our analysis are obtained by 
diagonalizing corresponding mass matrix numerically. The necessary amplitudes of tree 
and one-loop Feynman diagrams 2 and relevant effective couplings in the model are derived 
analytically by hand. By encoding the derived analytical formulas of the decay amplitudes 
into the numerical program, values of total amplitude and hence decay rate can be obtained. 
In the computation of the total decay width of light Higgs, we include all significant decay 
channel such as bb, t~t + , WW*, ZZ*, 77, and gg, as well as the RPV decay h° —> fif + rfi. 
With the RPV partial decay width rate for the channel and total decay width, the branching 
ratio can be obtained. 

During numerical analysis, following restrictions and assumptions are used: |/i |, \A U \, 



1 Higgses mix with sleptons via RPV terms, but we can still identify Higgses out of other sleptons due to 
the foreseeable smallness of RPV parameters. 

2 During numerical computation of Feynman diagrams, LoopTools package is used for the evaluation of loop 
integrals [2^. 
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| Ai| and | A x \ < 2500GeV; A e is set to be zero since its influence is quite small and negligible; 
rh 2 L (without zeroth component) < (2500GeV) 2 with off-diagonal elements zero; tan/3 

epton mass is demanded 
Q are pointing out a 



ranges from 3 to 60; heavy Higgs/sneutrino and charged Higgs/s 
to be above 200GeV. As discussed before, ATALS and CMS data 
Higgs-like boson with a mass around 125-126 GeV. After including the uncertainties on the 
experimental Higgs mass and loop corrections estimation to Higgs mass matrix in program, 
we allow the light Higgs mass to be in the range of 123 to 127GeV. For simplicity, we 
use the relation M 2 = = 2M\ between three gaugino masses and the condition that 

squarks of the first two families cannot be lighter than about O.8M3 . Therefore we take 
soft supersymmetry breaking scalar masses rhn = rhfj = m 2 D = (O.8M3 x identity matrix) 2 
for simplicity and M 2 < 2500GeV as boundary for analysis. The parameter setting is in 



24j . for instance. 



accordance with the gravity-mediated supersymmetry breaking picture 
We are interested in using a concrete setting that is compatible with known constraints but 
not otherwise too restrictive, to illustrate what we expect to be more generic features of the 
RPV signature. Besides, to prevent Bi from being unreasonably large, we set by hand 1% of 
Bq as upper bound in the circumstance of extraordinary loose or no available experimental 
bounds. 

We highlight below the most interesting combinations of RPV parameters giving signifi- 
cant contributions to the decays. 



A. Contributions from Bi fij combinations 

Among all Bifij combinations, only B 2 fi 3 an d B s fi 2 are expected to give important con- 
tribution to h° — > /if + tJjl because they contribute already at tree level (see diagram FIG. 
1 (Left Panel)) and B%\i 2 is enhanced by tau Yukawa coupling y e3 via a term proportional 
to y^MlB^y [(fioM 2 — sin2 / 9)M 2 ] thus becomes largest one among all B^/s (M s 
denotes a generic real scalar mass eigenvalue). It also happens to B 2 fi 3 , but with a muon 
Yukawa y e2 instead. These two combinations get their best values under small fiQ and M 2 
as can be seen from the formula. 

On the other hand, the values of Bi and Bi^ij are highly constrained separately by their 
loop contribution to neutrino mass matrix [25] ; a non-zero fij will induce tree level neutrino 
mass, hence it is constrained. In the meantime, leptonic radiative decays like fi — > etc. 
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also give upper bounds on Bifij, say, l-B^A^I and 1-^3/41 c2 10~ 4 I Mo 1 3 0I- 

By requiring tree level neutrino mass and all the loop contribution from Bi and Bifij 
to neutrino mass matrix elements within leV, we can get B.R. (-B2/X3) = 1 x 10~ 15 and 
B.R.(53yU2) = 1 x 10~ 13 at most. In this situation, restrictions by leptonic radiative decays 
are relatively loose such that we don't need to bother with it. 

If we allow neutrino masses being larger than leV while making /i large enough to 
unfasten restrictions by leptonic radiative decays, the main constraint will come from B^fij 
loop contribution. By demanding -82^3 and B^Hi contribution to the corresponding neutrino 
mass terms to be only smaller than 190keV, the possible branching ratio for each combination 
increases significantly to 9 x 10~ 6 and 7 x 10~ 4 , respectively. Note that B2 or B% alone can 
give contributions to h° — > \lt + tJX decay as well. Such contributions are included in the 
Bi/ij result (and -BjA, BiA x thereafter). 



B. Contributions from Bi A combinations 

Among all SjA's, it is obvious that -B1A123, -B 1 A 132 , -B2A232 and -B3A233 are the most 
important ones because they can contribute to the amplitudes at the tree level already. This 
tree level contribution can be approximated by .Mtree ~ -BjA(tan/3 sin a — cos a)/ [y2M^\ 
where a is the mixing angle between the two CP-even neutral Higgs. 



Table 1. Bi A contributions to branching ratio of h° — > fir + 77/ 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


Bi A123 


1 x 1(T 5 


4 x 10" 5 


Bi A 132 


3 x 


7 x 1(T 5 


B2 A232 


3 x 1(T 5 


6 x KT 2 


B3 A233 


3 x 


3 x KT 2 



The value of Bi is constrained by B^Bj neutrino mass loops [25(, while that of A by 
charged current experiments 26|]. Generally speaking, increasing supersymmetry breaking 
slepton mass and gaugino mass leads to heavier charged slepton, sneutrino and neutralino 
hence raises upper bounds for Bi and A . Due to the strong dependence upon Bi and A by 
tree level amplitude, it is still favorable even though heavy sneutrino tend to suppress the 
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branching ratio. Besides, /xq should not be too small in order to make the product of Bi and 
A be below the bounds from leptonic radiative decays, i.e. | -S* Ai 32 1 , |-E>iA^ 23 |, 1-02-^232 1 and 
I.B3A233I J$ 1.4 x 10~ 3 |/!o| 2 [12]. Even under stringent neutrino mass < leV constraint, we 
find that the four combinations could give sizeable branching ratios of about 10 -5 , which is 
large enough to be possibly probed at LHC experiments. 

If we do not impose the sub-eV neutrino mass bounds, B 1 X 12 3, B 1 X i32 and -B2A232 still 
get constrained by BiBj neutrino mass loops, while -B3A233 is limited by our manual "1% of 
B " bound. This causes their different magnitude of amplitudes. Comparing with neutrino 
mass < leV case, branching ratios from B\ raise a bit, while that from B 2 and B 3 increase 
significantly by 3 orders of magnitude and reach about 10~ 2 . 



C. Contributions from Bi A combinations 

The contributions may be quite interesting because this will be like the first experimental 
signature of the RPV A-parameters, and they are not constrained by the radiative decays 



12j. However, A x only plays its role in a single loop diagram (FIG. 3 (Left Panel)) via neutral 
scalar-charged scalar-charged scalar (h°<f) + <f)~) coupling. It is expected to give important 
contributions for low charged scalar mass to avoid strong suppression. A x do not have 
experimental constraints 3 , and can consequently take any value. But Bi is still limited by 
neutrino mass loops as before. Similar to BiX case, when considering experimental neutrino 
mass, different mass constraints induce different magnitude of branching ratios (for B\A X 
and B 2 A X ) while B 3 is constrained by the "1% of Bq" bound. 

Table 2. BiA x contributions to branching ratio of h° — > /jt + r/2 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


Bi A x 23 


5 x KT 11 


2 x KT 10 


Bi A 132 


5 x KT 11 


2 x KT 10 


B 2 A 232 


5 x KT 11 


7 x KT 7 


B3 A x 33 


5 x KT 11 


1 x KT 7 



3 Recall that under the parameterization adopted here, the Li superfields all have zero VEV. Hence, A x as 
defined does not contribution, for example, to b — » 57 at one-loop level. 
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In our setting, branching ratios from B^A X can at most reach the order of 1CT 7 . However, 
if we allow A x to be as large as 3.1TeV (or even higher), branching ratio of the order of 1CT 6 
and above is possible. Since decay rate is proportional to amplitude square and hence A x 
square, it is easy to see how branching ratio modifies as A x increases as can be seen in Fig. 
3 (Right Panel). 

IV. SUMMARY 

There are different scenarios to achieve LFV Higgs decay in supersymmetric standard 
model depending on how neutrino masses are implemented. We have shown that even 
with RPV parameters only, significant contributions to h° — >■ /if + r/2 are possible. In our 
analysis, Bi A's undoubtedly provide largest branching ratios. Even with stringent neutrino 
mass constraint < leV, several combinations of Bi A can still give branching ratios beyond 
1CT 5 which should not be overlooked in future collider experiments. We also listed branching 
ratios as references in case that the sub-eV neutrino mass is relaxed. Such a case is unlikely 
but not definitely ruled out. 

At the end, we would like to encourage our experimentalist collegues at CMS and ATLAS 
to investigate those flavor violating Higgs decays because they may give some complementary 
informations about lepton flavor violating couplings, which are also relevant for neutrinos 
physics. Moreover, a typical cross section of MSSM 125 GeV at 8 TeV energy is of the order 
lOpb. For Br(h° — >■ /if + rfi) of the order 10" 4 and with a luminosity of the order 10 fb -1 , 
this would lead to 10 raw /if + tJjl events with almost no SM background. This would get 
amplified with the 14 TeV energy for future LHC run which may give an even better chance 
to probe the lepton flavor violating couplings. 
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FIG. 1 (Left Panel): An example of tree diagrams. Charged lepton \i transforms into 
higgsino via RPV parameter fi2, while light Higgs transforms into sneutrino z>3 via B%. 
This diagram gives major contribution to h° — > tJl. 

(Right Panel): Branching ratio from -B3/Z2 as a function of B 3 and H2, with M 2 = 850GeV, 
ml u = m 2 Eii = (800GeV) 2 , /i = 2500GeV = A u = -A d , tan/3 = 60. m A 200 to 320GeV 
in the permitted region. Solid red line [m v bound) comes from demanding 23-element of 
neutrino mass matrix <190keV, while solid vertical line indicates bound from Higgs mass. 
If Bj, increases across the bound line, the mass of Higgs (both light and heavy ones) will 
be below permitted lower bounds. 



14 



1.5 



0.5 





1 -1 1 

%232 bound 




- 


• 

V 

Br=10 


- 




B 2 bound 






10 6 




- 


10 7 






10" 8 






10 9 - 






I — 1 1— 





500 



1000 
B 2 (GeV) 2 



1500 



2000 




1000 2000 3000 4000 5000 6000 7000 8000 

■\2 



B 2 (GeV)* 



FIG. 2 (Left Panel): Branching ratio from B 2 X 232 , with M 2 = 2500GeV, m 2 Lii = rh% u = 
(2500GeV) 2 , /i = 1900GeV = A u = -A d , tan/3 = 60. M A = 200 to 205GeV in the permit- 
ted region. Upper bounds of B 2 and A 232 significantly raised by large soft breaking masses 
rh 2 Lii , rh 2 Eii and gaugino mass M 2 . Solid vertical line (B 2 bound) comes from demanding 
22-element of neutrino mass matrix <leV. 

(Right Panel): Branching ratio from B 2 A232 without sub-eV neutrino mass constraint, with 
M 2 = 850GeV, m 2 Lii = m 2 m = (1400GeV) 2 , /i = 2500GeV = A u = -A d , tan/3 = 60. 
Ma — 200 to 390GeV in the permitted region. Solid vertical line (B 2 bound) comes from 
demanding 22-element of neutrino mass matrix <190keV. 
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FIG. 3 (Left Panel): The only one-loop diagram in which A x has its stage. A x participates 
in via neutral scalar-charged scalar-charged scalar coupling. 

(Right Panel): Branching ratio from B2A X 32 as a function of B2 and ^232' with M2 = 
850GeV, m 2 Lii = m 2 Eii = (555GeV) 2 , /i = 2500GeV = A u = -A d , tan/3 = 60. M A ^ 200 
to 300GeV in the permitted region. Solid vertical line comes from demanding 22-element of 
neutrino mass matrix <190keV. 



